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We study the effect of color superconductivity and neutrino trapping on the deconfinement transi- 
tion of hadronic matter into quark matter in a protoneutron star. To describe the strongly interacting 
matter a two-phase picture is adopted. For the hadronic phase we use different parameterizations 
of a non-linear Walecka model which includes the whole baryon octet. For the quark matter phase 
we use an SU (3) / Nambu-Jona-Lasinio effective model which includes color superconductivity. We 
impose color and flavor conservation during the transition in such a way that just deconfined quark 
matter is transitorily out of equilibrium with respect to weak interactions. We find that deconfine- 
ment is more difficult for small neutrino content and it is easier for lower temperatures although 
these effects are not too large. In addition they will tend to cancel each other as the protoneutron 
star cools and deleptonizes, resulting a transition density that is roughly constant along the evo- 
lution of the protoneutron star. According to these results the deconfinement transition is favored 
after substantial cooling and contraction of the protoneutron star. 

PACS numbers: 12.39.Fe, 25.75.Nq, 26.60.Kp 



I. INTRODUCTION 



It is currently a matter of speculation the actual occurrence of quark matter during protoneutron star (PNS) 
evolution. The standard scenario for the birth of neutron stars indicates that these objects are formed as consequence 
of the gravitational collapse and supernova explosion of a massive star [Ml. Initially, PNSs are very hot and lepton- 
rich objects, where neutrinos are temporarily trapped. During the first tens of seconds of evolution the PNS evolves 
to form a cold (T < 10 10 K) catalyzed neutron star [H-Q- As neutrinos are radiated, the lepton - per - baryon content 
of matter goes down and the neutrino chemical potential tends to essentially zero in ~ 50 seconds [1[ . Deleptonization 
is fundamental for quark matter formation inside neutron stars, since it has been shown that the presence of trapped 
neutrinos in hadronic matter strongly disfavors the deconfinement transition 0, @- In fact, neutrino trapping makes 
the density for the deconfinement transition to be higher than in the case of neutrino-free hadronic matter. As 
a consequence, the transition could be delayed several seconds after the bounce of the stellar core. However, the 
calculations presented in 0, [El were performed employing the MIT Bag model for the description of quark matter and 
did not include the effect of color superconductivity. As we shall see in the present work, the use of the Nambu-Jona- 
Lasinio model and the inclusion of color superconductivity may change qualitatively the effect of neutrino trapping 
in the deconfinement conditions. 

As emphasized in earlier works [il-fllj|. an important characteristic of the deconfinement transition in neutron stars, 
is that just deconfined quark matter is transitorily out of equilibrium with respect to weak interactions. In fact, 
depending on the temperature, the transition should begin with the quantum or thermal nucleation of a small quark- 
matter drop near the center of the star. On the other hand, the flavor composition of hadronic matter in /^-equilibrium 
is different from that of a /3-stable quark-matter drop. Roughly speaking, the direct formation of a /3-stable quark-drop 
with N quarks will need the almost simultaneous conversion of ~ A/3 up and down quarks into strange quarks, a 
process which is strongly suppressed with respect to the formation of a non /3-stable drop by a factor ~ G Fcrmi . For 
typical values of the critical-size /3-stable drop (N ~ 100 — 1000 joT) the suppression factor is actually tiny. Thus, 
quark flavor must be conserved during the deconfinement transition [4-llOj . When color superconductivity is included 
together with flavor conservation, the most likely configuration of the just deconfined phase is 2SC provided the 
pairing gap is large enough [9|. The relevance of this 2SC intermediate phase (a kind of activation barrier) has been 
analyzed for deleptonized neutron stars [Io|, [HI but not for hot and lepton-rich objects like PNSs. 

In the present paper we shall analyze the deconfinement transition in protoneutron star conditions employing the 
Nambu-Jona-Lasinio model in the description of quark matter. For the hadronic phase we shall use a model based 
on a relativistic Lagrangian of hadrons interacting via the exchange of <r, p, and u mesons [l3l |. For simplicity, the 
analysis will be made in bulk, i.e. without taking into account the energy cost due to finite size effects in creating a 
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TABLE I: Parameters of the hadronic equation of state. For each parametrization we give the maximum mass M max of a 
hadronic star. 



drop of deconfined quark matter in the hadronic environment. 

The present article is organized as follows. In Sec. II we present the main aspects of the non-linear Walecka model 
describing the hadronic phase. In Sec. Ill we present the generalities of the model we use for the quark phase. In 
Sec. IV we show our numerical results and finally in Sec. V we discuss our results and present the conclusions. 



II. THE HADRONIC PHASE 

For the hadronic phase we shall use a non-linear Walecka model (NLWM) [l3l - fil| which includes the whole baryon 
octet, electrons and electron neutrinos in equilibrium under weak interactions. The Lagrangian of the model is given 
by 

£ = Cb + Cm + C l , (1) 
where the indices B, M and L refer to baryons, mesons and leptons respectively. For the baryons we have 
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with B — n, p, A, £ + , £ , £ , S , and S°. The contribution of the mesons er, w and p is given by 
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where the coupling constants are g G B — %aB g<j, g^B — %uB g^ and g P B = x p b g P - The ratios x a s, i u b and x p b are 
equal to 1 for the nucleons and acquire different values for the other baryons depending on the parametrization (see 
Table I) . The leptonic sector is included as a free Fermi gas of electrons and electron neutrinos in chemical equilibrium 
with the other particles. 

There are five constants in the model that are determined by the properties of nuclear matter, three that determine 
the nucleon couplings to the scalar, vector and vector-isovector mesons gol^o, 3u/to w , g p /m p , and two that determine 
the scalar self interactions b and c. It is assumed that all hyperons in the octet have the same coupling than the 
A. These couplings are expressed as a ratio to the nucleon couplings mentioned above, that we thus simply denote 
x a , x u and x p . In the present work we use two parameterizations for the constants. One of them is the standard 
parameterization GM1 given by Glendenning-Moszkowski [13j . as shown in Table I. This parametrization employs 
"low" values for x a , x u and x p . The parametrization GM4 employs larger values of these couplings. This makes the 
EOS stiffer and increases the maximum mass of hadronic stars to 2.2 Mq, see Table I. 

The derivation of the equations describing the model is given in detail in (l6| . The total pressure P and mass - 
energy density p are given by: 

2 1 / \ -2 -, -, t / \ 2 



Here P$ and p$ are the expressions for a Fermi gas of relativistic, non-interacting particles: 

* = \wf I d3p wrifw {MT) + mi (6) 
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Pi = WrJ d3p {p2 + mf)1/2 iMT) + /i(T)) ' (7) 

where fi(T) and fi(T) are the Fermi - Dirac distribution functions for particles and antiparticles respectively: 

fi(T) = (exp([{p 2 + to* 2 ) 1/2 tf]/T) + l)" 1 , (8) 

MT) = (exp([(p 2 + mff/ 2 + »*}/T) + l)" 1 . (9) 

Note that for baryons we use, instead of masses To; and chemical potentials fa, "effective" masses to* and chemical 
potentials p* given by: 

to* — Tfii -\- x al (g a a), (10) 

A*< = Mi - ^ Pb - V 7 3i (j^) P/ 3 > ( n ) 

where I^i is the third component of the isospin of each baryon. 

The weighted isospin density p l3 and the weighted baryon density p B are given by: 

Pl 3 = ^2 X Pi I 3i n ii ( 12 ) 
i=B 

Pb = ^2 x ^^> ( 13 ) 

i=B 

being rn the particle number density of each baryon: 



n, = 



(2tt) 3 

The mean field g a o satisfies the equation: 



J cPpifiW-MT)). (14) 
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(g a a) + bm n {g a a) 2 + c(5 CT cr) 3 = ^ x oi n \^ ( 15 ) 

i=B 

where n\ is the scalar density: 

nt = Wfl d3p WT<W {MT) + MT)) - (16) 

The hadron phase is assumed to be charge neutral and in chemical equilibrium under weak interactions. Electric 
charge neutrality states: 

n P + — n s- — n E- — n e = 0. (17) 

Chemical weak equilibrium in the presence of trapped electron neutrinos implies that the chemical potential fa of 
each baryon in the hadron phase is given by: 

Pi = qsfai - q e (pe - p Ve ), ( 18 ) 

where qs is its baryon charge and q e is its electric charge. For simplicity we are assuming that muon and tau neutrinos 
are not present in the system, and their chemical potentials are set to zero. 

All the above equations can be solved numerically by specifying three thermodynamic quantities, e.g. the temper- 
ature T, the mass-energy density p and the chemical potential of electron neutrinos in the hadronic phase p^ . 
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III. THE QUARK MATTER PHASE 

In order to study the just deconfined quark matter phase we use an SU (3) / NJL effective model which also includes 
color superconducting quark-quark interactions. The corresponding Lagrangian is given by 



C 
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where rh = Aia,g{m u ,m c i^m s ) is the current mass matrix in flavor space. In what follows we will work in the isospin 
symmetric limit m u = wid = m. Moreover, r, and Aj with i = 1, .., 8 are the Gell-Mann matrices corresponding to the 
flavor and color groups respectively, and To = -\/2/3 1/. Finally, the charge conjugate spinors are defined as follows: 
ipc = C ip T and ipc = ^ T C, where ip = ip^-f is the Dirac conjugate spinor and C = ij 2 j°. 

To be able to determine the relevant thermodynamical quantities we have to obtain the grand canonical thermo- 
dynamical potential at finite temperature T and chemical potentials /i/ c . Here, / = (u, d, s) and c = (r, g, b) denotes 
flavor and color indices respectively. For this purpose, starting from Eq. (|19[) . we perform the usual bosonization of 
the theory. This can be done by introducing scalar and pseudoscalar meson fields a a and 7r a respectively, together with 
the bosonic diquark field A a- In this work we consider the quantities obtained within the mean field approximation 
(MFA). Thus, we only keep the non- vanishing vacuum expectation values of these fields and drop the corresponding 
fluctuations. For the meson fields this implies a = a a T a = diag(a„, ad, a s ) and 7r Q = 0. Concerning the diquark mean 
field, we will assume that in the density region of interest only the 2SC phase might be relevant. Thus, we adopt the 
ansatz A5 = A7 = 0, A2 = A. Integrating out the quark fields and working in the framework of the Matsubara and 
Nambu-Gorkov formalism we obtain the following MFA quark thermodynamical potential (a detailed procedure of 
calculation can be found in Refs. [17H19l | ) 
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where A is the cut-off of the model and uj(x, y) is defined by 
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Vl = Hub , V2 = fJ-db , V3 = fJ-sr , 2/4 = fJ>sg , 2/5 = Vsb , 2/6,7 = ^ , 2/8,9 = » • (22) 



Here, E = \Jk 2 + M 2 and E s = yjk 2 + M 2 , where Mf — mj + 07. Note that in the isospin limit we are working 
c« = &d — a an d, thus, M u = Md = M. 

The total thermodynamical potential of the quark matter phase (QMP) is obtained by adding to Qmfa the 
contribution of the leptons. Namely, 

n QMP (T,fi fc , fi e , n„ e ,a,a s , I A|) = Q.f FA (T, » fc , a, a s , | A|) + O e (T, /i e ) + (l v . (T, Ml/e ) - f> vac (23) 

where Q e and Cl„ e are the thermodynamical potentials of the electrons and neutrinos, respectively. For them we use 
the expression corresponding to a free gas of ultra-relativistic fermions 



\ 24tt 2 12 360 



(24) 



where I = e, v e and the degeneracy factor is 7 e = 2 for electrons and 7^ = 1 for neutrinos. Notice that in Eq. (|2"3")l 
we have subtracted the constant Sl vac in order to have a vanishing pressure at vanishing temperature and chemical 
potentials. 

From the grand thermodynamic potential Qqmp we can readily obtain the pressure P — —CIqmp, the number 
density of quarks of each flavor and color 7if c = —9ClQMp/9fif c , the number density of electrons n e = —dQQMp/d/-i e , 
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TABLE II: The two sets of NJL parameters. 



and the number density of electron neutrinos n„ e — —dQQMp/d(J, Ue . The corresponding number densities of each 
flavor, nf, and of each color, n c , in the quark phase are given by nf — X) c n /c an< ^ n c ~ S/ n /c respectively. The 
baryon number density reads tib = | ^2f c n fc = ( n u + «d + n s )/3. Finally, the Gibbs free energy per baryon is 



1 



ffquark = ( ^ Vfc ™/c + Me ™e + V-v* n Ve | . (25) 

fc 



For the NJL model we use two sets of constants shown in Table HH The sets 1 and 2 were taken from [2(| and [2l[ 
respectively, but without the 't Hooft flavor mixing interaction. The procedure, obtained from [22j is to keep A and 
m fixed, then tune the remaining parameters G and m s in order to reproduce M — 367.6 MeV and M s = 549.5 MeV 
at zero temperature and density. The resulting parameter sets are given in Table HU 

In order to derive a quark matter EOS from the above formulae it is necessary to impose a suitable number of 
conditions on the variables {/i/ c }, £te, Hv e , &, &s and A. Three of these conditions are consequences from the fact that 
the thermodynamically consistent solutions correspond to the stationary points of fl with respect to a, <j s , and A. 
Thus, we have 

dn QM p/da = , dn QMP /da s =a , dCl QMP /d\A\ = 0. (26) 

To obtain the remaining conditions one must specify the physical situation in which one is interested in. As in 
previous works [34Toj . we are dealing here with just deconfined quark matter that is temporarily out of chemical 
equilibrium under weak interactions. The appropriate condition in this case is flavor conservation between hadronic 
and deconfined quark matter. This can be written as 

Yf =yf f = u,d,s,e,v e (27) 

being Yp = / rig and = nj / the abundances of each particle in the hadron and quark phase respectively. 
In other words, the just deconfined quark phase must have the same "flavor" composition than the /3-stable hadronic 
phase from which it has been originated. Notice that, since the hadronic phase is assumed to be electrically neutral, 
flavor conservation ensures automatically the charge neutrality of the just deconfined quark phase. The conditions 
given in Eq. (|27p can be combined to obtain 

rid = £ n u , n s = i] n u , n Ve — n n u , 3n e = 2n u - n d - n s , (28) 

where m is the particle number density of the z-species in the quark phase. The quantities £ = Yj 1 /Y^ , -q = Y^ /Y^ 
and k = Yj 1 /Y^ are functions of the pressure and temperature, and they characterize the composition of the hadronic 
phase. These expressions are valid for any hadronic EOS. For hadronic matter containing n, p, A, S + , S , S~ 
and 5°, we have 

^ _ n p + 2n n + n A + n s o + 2n s - + n s - 
2n p + n n + nA + 2n s + + n s o + n s o ' 

n\ + n s + + n s o + n s - + 2n H o + 2n s - 

V = — ^ ^ 1 (30) 

2n p + n n + n A + 2n s + + n s o + n s o 

(31) 



up 



2n p + n n + n\ + 2n s + + n s o + n s o 



Additionally, the deconfined phase must be locally colorless; thus it must be composed by an equal number of red, 
green and blue quarks 



n r = n g = rib- 



(32) 
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Also, ur, ug, dr, and dg pairing will happen provided that |A| is nonzero, leading to 

™ur — Tldg 3 rh u g — Tldr- (^3) 

In order to have all Fermi levels at the same value, we consider Q 

Ti U g Tl ur , f^sb Tlsr- (^^) 

These two equations, together with Eqs. (1321) and (133[) imply that n ur = n ug = n^r = n^g and n sr = n sg = n s i, Q. 

Finally, including the conditions Eqs. (|26j) we have 13 equations involving the 14 unknowns (a, a s , |A|, p e , p Ue and 
{uf c }). For given value of one of the chemical potentials (e.g. p ur ), the set of equations can be solved once the values 
of the parameters £, ry, k and the temperature T are given. Instead of p ur , we can provide a value of the Gibbs free 
energy per baryon £fq U ark or the pressure P and solve simultaneously Eqs. (|28l) - (p4|) together with Eq. (|26|) in order 
to obtain tr, cr s , |A|, p e , p U(s and {/i/ c }. 



IV. DECONFINEMENT TRANSITION IN PROTO-NEUTRON STAR MATTER 



In order to determine the transition conditions, we apply the Gibbs criteria, i.e. we assume that deconfinement will 
occur when the pressure and Gibbs energy per baryon are the same for both hadronic matter and quark matter at a 
given common temperature. Thus, we have 

g H = g Q , P H = P Q , T H = , (35) 

where the index H refers to hadron matter and the index Q to quark matter. According to these conditions (together 
with the equations of Sections 2 and 3), for a given temperature T H and neutrino chemical potential of the trapped 
neutrinos in the hadronic phase , there is an unique pressure P at which the deconfinement is possible. Instead of 
P, we may characterize the transition point by giving the Gibbs free energy per baryon g, or alternatively, the mass- 
energy density of the hadronic phase pn (see Figs. 1-3). We emphasize that, according to the present description, P 
and g are the same in both the hadronic phase and the just deconfincd phase. However, the mass-energy density pn 
and pq at the transition point are different in general. Similarly, while the abundance Y Vg of neutrinos is the same in 
both the hadronic and just deconfined quark phases, the chemical potentials p^ and p^ are different. 

According to numerical simulations [lH3j], during the first tens of seconds of evolution the protoneutron star cools 
from T ~ 40 MeV to temperatures below 2-4 MeV. In the same period, the chemical potential p^ of the trapped 
neutrinos evolves from ~ 200 MeV to essentially zero. Thus, in order to consider typical PNS conditions we have 
solved Eqs. (|35|) together with the equations of Sections 2 and 3 for temperatures in the range — 60 MeV and p^j in 
the range — 200 MeV. The results are displayed in Figs. 1-3 for all the parameterizations of the equations of state 
given in previous sections. 

In Fig. [T]we display the results for the GM1 parametrization of the hadronic EOS. In the left panel of Fig. [T]we show 
the Gibbs free energy density per baryon g at which deconfinement occurs versus the temperature T for three different 
values of the neutrino chemical potential in the hadronic phase [p^ e = 0, 100, 200 MeV). In the right panel the same 
results are shown but as a function of the mass-energy density of the hadronic phase (in units of the nuclear saturation 
density po = 2.7 x 10 14 g cm 3 ). In both figures, if the thermodynamic state of hadronic matter (characterized e.g. 
by T H , ph and p^J ) lies to the left of the curve corresponding to the same p^ , then the deconfinement transition is 
not possible. In the right side region of a given curve the preferred phase is deconfined quark matter. Notice that the 
transition's Gibbs free energy is an increasing function of p^ . However, the transition density of the hadronic phase 
slightly decreases as p^j increases. In Fig. [2] we display the results for the GM4 parametrization of the hadronic EOS. 
The results are qualitatively the same but the transition densities are smaller than those for GM1 by ~ 30%. 

In Fig. 3 we show the behavior of the transition's density as a function of the chemical potential of trapped neutrinos 
p^ for two specific temperatures (T = 2 and 30 MeV). It is clearly seen that for a fixed temperature the effect of 
deleptonization is to inhibit the transition. This effect is not very large; at fixed temperature there is a slight increase 
by less than a 10% when p^ falls from 200 MeV to MeV. On the other hand, the effect of cooling works in the 
opposite direction because pairing tends to help the transition and the gap increases as the temperature goes down. 
The effect of cooling is also small; at fixed p^ there is a slight decrease of the transition density pn by less than a 
10% when the temperature falls from 30 MeV to 2 MeV. Both effects tend to cancel each other as the PNS cools and 
deleptonizes, resulting a transition density that is roughly constant along the evolution of the protoneutron star. 
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FIG. 1: Left panel: the Gibbs free energy density per baryon g at which deconfinement occurs versus the temperature T for 
three different values of the neutrino chemical potential in the hadronic phase (fi,„ e = MeV in dashed line, /x^_ = 100 MeV in 
dotted line and /x^, = 200 MeV in full line). Right panel: the mass-energy density of the hadronic phase at which deconfinement 
occurs versus the temperature T, for the same values of /x* given in the left panel (density is given in units of the nuclear 
saturation density po )• The hadronic phase is described by the GM1 parametrization of the EOS. For the quark phase we 
adopt the two parameterizations of the NJL model given in Table II. In both figures, if the thermodynamic state of hadronic 
matter (characterized by { T H , g 11 , /U^,} or by {T H , pn, t 1 ^ } ) lies to the left of the curve corresponding to the same //^, 
then the deconfinement transition is not possible. In the right side region of a given curve the preferred phase is deconfined 
quark matter. 




FIG. 2: Same as Fig. 1 but employing the GM4 parametrization of the hadronic equation of state. 



V. CONCLUSIONS 



In this paper we have investigated the role of color superconductivity in the deconfinement transition from hadronic 
matter to quark matter at finite temperature and in the presence of a trapped neutrino gas. The study presented 
here is relevant for the first tens of seconds of evolution of newly born protoneutron stars. 

In our analysis we used a two phase description where we employed the Nambu-Jona-Lasinio model in the descrip- 
tion of quark matter (Sec. Ill) and a non-linear Walecka model which includes the whole baryon octet, electrons 
and electron neutrinos in equilibrium under weak interactions in the description of hadronic matter (Sec. II). Decon- 
finement is assumed to be a first order phase transition and the just deconfined quark phase is assumed to have the 
same "flavor composition" than the /3-stable hadronic phase from which it has been originated (see [l2j and references 
therein). When color superconductivity is included together with flavor conservation [9|, the most likely configuration 
of the just deconfined phase is 2SC provided the pairing gap is large enough. This just deconfined phase is out of 
chemical equilibrium under weak interactions and thus it is very short lived but it is a a kind of "activation barrier" 
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FIG. 3: The mass-energy density of the hadronic phase at which deconfinement occurs as a function of the chemical potential 
of trapped neutrinos p^ . Results are given for two temperatures: T — 2 MeV in full line and T = 30 MeV in dashed line. We 
employed the GM1 (left panel) and the GM4 (right panel) parametrization of the hadronic EOS. Notice that there is a small 
decrease of the transition density pu for large fjj^ e . 



that determines the onset of the deconfinement transition. 

The main result of the present paper is that, within the NJL model, deconfinement is more difficult for small 
neutrino content and it is easier for lower temperatures. This effect is not very large, at least for the here-used 
parameterizations of the NJL model. At fixed temperature there is a slight increase by less than a 10% when pJJ 
falls from 200 MeV to MeV (see Fig. 3). The effect of cooling is also small; at fixed /i^ there is a slight decrease of 
the transition density pn by less than a 10% when the temperature falls from 30 MeV to 2 MeV (see Fig. 3). This 
is due to the fact that the pairing gap becomes larger as the temperature decreases and therefore the increase of the 
condensation term favors the transition at low temperatures. Both effects tend to cancel each other as the PNS cools 
and deleptonizes, resulting a transition density that is roughly constant along the evolution of the protoneutron star. 

The here-found behavior is qualitatively opposite to what it was found within the MIT bag model. In fact, previous 
analysis without including the effect of color superconductivity 0, [B| show that the presence trapped neutrinos pushes 
up the transition density to values much larger than for neutrino free matter. It was also found in [J, |5| that the 
transition is easier for larger temperatures. More recent results including the effect of color superconductivity within 
the MIT Bag model [23| show that the transition density increases with neutrino trapping but (in coincidence with 
the here found results) the pairing gap favors the transition as the temperature decreases. 

In spite of some differences between the results within the NJL and the MIT bag model description of quark 
matter some general conclusions may be obtained about the effect of color superconductivity in the deconfinement 
transition. First, when color superconductivity is present the deconfinement density is not so strongly affected by 
neutrino trapping as it is in the unpaired case. Second, color superconductivity makes the transition easier at lower 
temperatures and the dependence of the deconfinement density with T is much smaller than in the unpaired case. 

During cooling and deleptonization of the protoneutron star the temperature and the chemical potential of trapped 
neutrinos fall abruptly in a few seconds and there is also some contraction of the whole neutron star. It is interesting 
to note that although the density increase is not too large, it may be comparatively important for the deconfinement 
transition because the effects of temperature and neutrino trapping are smoothed by color superconductivity. Accord- 
ing to our results the deconfinement transition is favored after substantial cooling and contraction of the protoneutron 
star but full numerical simulations of protoneutron star evolution are needed in order to determine whether and when 
the deconfinement conditions are attained. 
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